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Abstract 

We propose the tube model as a first step in solving the bound 
state problem in light-front QCD. In this approach we neglect trans- 
verse variations of the fields, producing a model with 1+1 dimensional 
dynamics. We then solve the two, three, and four particle sectors of 
the model for the case of pure glue SU(3). We study convergence to 
the continuum limit and various properties of the spectrum. 

The ultimate goal of light-front field theory is to start with the QCD 
Lagrangian and, with a minimum of approximation, calculate the hadron 
spectrum. The basic idea behind this approach is to use Hamiltonian tech- 
niques in the coordinate system best suited for relativistic dynamics H . We 
believe that the "tube model" provides a good starting point towards that 
goal. Simply put, we take full 3+1 dimensional QCD and restrict our atten- 
tion to a region that is small in the two transverse directions and is large in 
the longitudinal dimension, hence the name "tube model." Since the trans- 
verse region is small, we neglect any transverse derivatives of the fields as 
an approximation. The resulting theory is a 1+1 dimensional gauge theory 
coupled to adjoint matter. We then employ light-front Hamiltonian tech- 
niques to produce the low energy spectrum of the theory. This approach 
was first explored in the transverse lattice paper of Bardeen, Pearson, and 
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Rabinovici [0] (that is, one can think of the tube model as two links on a 
transverse lattice) and has become increasingly popular in the last few years 
under the name "dimensional reduction" [| . 

The tube model has the virtues of a 1+1 dimensional theory, in that it 
is superrenormalizable and computationally tractable, while retaining some 
of the dynamics of the full theory. Thus, as we demonstrate in this paper, 
one can hope to understand the behavior of one tube quite thoroughly. The 
subsequent step will be much more difficult: to couple several of these tubes 
together to make a full 3+1 dimensional theory. 

In this paper, we study only the pure glue sector of the theory. Of course, 
an important next step will be to include fermions. Although one may impose 
antiperiodic boundary conditions on the fermion fields, a consistent formula- 
tion that includes fermions must impose periodic boundary conditions on the 
gauge fields. In anticipation of this situation, we impose periodic boundary 
conditions for our gauge fields; this will have some important consequences 
for our numerical results. 

The paper is organized as follows. In Sec. |1] we discuss (and dismiss) the 
gauge zero mode. We define the tube model for pure glue in Sec. §. Then in 
Sees. |, £|, and || we discuss the two, three, and four particle sectors of the 
theory. In addition, our technique for measuring parity is demonstrated in 
Sec. f|. Finally, we discuss some aspects of our results in Sec. |6|. 

1 The gauge zero mode 

Before we begin, we must discuss the "gauge zero mode," V = J dx~ A_. 
As has been noticed by several authors in the context of both light-front || 
and equal time || quantization, one cannot set V = by choice of gauge 
because it is a genuine dynamical degree of freedom. Only the gauge choice 
d-A_ = is allowed. 

What are the effects of V on the tube model particle spectrum? In the 
Lagrangian, we must replace the longitudinal derivative cL by the covariant 
version cL + igV. If we regulate the theory by imposing periodic boundary 
conditions (DLCQ), the longitudinal momentum becomes discrete and we 
find that the effect of V is to shift momenta by an amount that is less than 
one longitudinal momentum spacing. Clearly, any direct effect of this mode 
on ordinary dynamical modes is negligible in the large volume limit L — > 
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oo. However, the zero mode sector of the theory may become qualitatively 
different. 

For instance, some of the zero modes that were otherwise constrained 
become dynamical. In light-front field theory, equations of motion are typi- 
cally of the form <9 + (<9_ + igV)(j) = U ((/>), where U(4>) is the potential term. 
When V — it is clear why zero modes <9_0 = are not dynamical: they 
have no conjugate momentum. However, note that when V is nonzero, the 
associated kinetic term is small (since V is small) relative to the kinetic term 
of any non-zero mode. Thus, one may just as well treat these zero modes as 
non-dynamical (that solution of a constraint equation). 

Without the gauge zero mode, Gaufi' law requires that all N 2 — 1 gener- 
ators Q a of the gauge group (defined in Eqn. (|?p) annihilate physical states. 
If V is included, only the N — 1 diagonal generators must annihilate physi- 
cal states. Consequently, the spectrum contains many more "color singlet" 
states if V is nonzero. However, one finds that these additional states have 
eigenvalues that diverge as one takes the continuum limit L — > oo and they 
do not affect the low energy spectrum [J], |j. 

Although V may be quite important for understanding the vacuum state 
itself, its effects on the particle spectrum seem to be small in the type of 
model we are studying. Thus, we set V = as a dynamical approximation. 
It is not clear that A_ can be so easily neglected in the full 3+1 dimensional 
theory. For instance, in 3+1 dimensions, A_ has a "proper zero mode," 
/ dx~ A-, that is constrained along with a dynamical "global zero mode," 
/ dx~ d 2 x r A_. The effects of the proper zero mode on the particle spectrum 
could be quite important [|]. 

2 The Hamiltonian 

In this section, we define the tube model for 3+1 dimensional SU(N) gauge 
theory, introduce our notation, and write out the light-front Hamiltonian. We 
introduce light front coordinates x^ = x ± = (x° ±a; 3 )/v / 2, where x + is light- 
front "time" and x~ is the longitudinal spatial coordinate; transverse indices 
are labeled by r,s G {1,2}. Next, we put the system in a box of dimen- 
sions (L, L±, L±) in the (x~ , x , x 2 ) directions and impose periodic boundary 
conditions, or "discretized light-cone quantization" (DLCQ) |], |10| . As our 
central assumption, we neglect variations of the fields in the transverse direc- 
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tions and set d r A^ = 0. Since this approximation violates transverse gauge 
invariance and rotational invariance, we are free to introduce a transverse 
mass term tr A r A r . As we shall see, inclusion of this term will be necessary 
for proper renormalization of the model. Thus the Lagrangian is 

£ g augc = L\ J dx~ tr + 2F +r F_ r - ^F rs F rs - /i 2 A r A r \ , (1) 
where 

F-+ = d^A + (2) 

F_ r = d^A r (3) 

F +r = d + A r + ig [A+, A r ] (4) 

F rs = ig [A r ,A s ] . (5) 

The A + field is a constrained degree of freedom. It obeys the constraint 
equation 

- ^<9_<9_,4 + = J+ = i [d^A r , A r ] , (6) 

where J 11 is the conserved current. The conserved charge associated with J 71 
is 

Q = L\J dx~ J + = %L\ J dx~ [d^A r , A r ) . (7) 

According to Eqn. (H), Q must vanish. In the usual manner, we demand that 
Q annihilates physical states. 

The generators of SU(N) obey the relations T a ,T b = if abc T c and 

tr T a T b = 5 ab /2. Writing A r = A a r T a ; the momentum conjugate to A a r is 
d-A%. Thus the zero mode of A a r has no conjugate momentum and is not 
a true dynamical degree of freedom. For the rest of this discussion, we will 
drop the zero mode associated with A r . Next, we construct the helicity op- 
erator J z which generates rotations in the (x x ,x 2 ) plane, the longitudinal 
momentum operator P + , and the Hamiltonian P~: 

J z = -AL 2 ± fdx-tTAid^A 2 , (8) 

P + = 2L\ j dx'tid^A r d.A r , (9) 



P- = Llj dx-tr^-g 2 J + -^yJ + + ^A r A r -^[A ri A s } 2 ^ . (10) 
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The inverse derivative operator is well-defined in this case since J + has no 
zero mode in the space of physical states. 

Now we are ready to quantize the theory. We can expand A r in terms of 
creation and annihilation operators: 



E E4r( 

A=±l fc=l V K 



a —X—2irikx / L < a f X 2irikx /L 



a k,\ € r e 



+ a k,X e r e 



The creation/annihilation operators obey the usual commutation relations, 



a k,Xi a LX> 



a fc,A' a LX' 



o o k> tdx,x> 

a k,X j a Z,A' 



. 



(12) 
(13) 



We define the polarization tensor e^: ef 1 = l/y/2 and ef 1 = ±i/y/2. Using 
the creation/annihilation operators and normal ordering, we find 



a fa 



a k,X a k,X 



k=l 



2ty 



P+ = Y.ka a k 



a f a 



Q a = -* E f ab ° a k,X a k,X 
k=l 



(14) 
(15) 

(16) 



Let us examine the Hamiltonian P~ in more detail. From the first two 
terms in Eqn. (|10D we obtain the four instantaneous interactions, the "self- 
induced inertias," and the mass term, 



P-> 



9 2 L 



E 



fabc fa'b'c 
m+n,m'+n' J J 



n' :nv :n,rn 



16n 2 L 2 , ^ 
9 2 L 



E 



\/mm'{m + n) 2 
■n a V V n b ' 

a n,-X a m,X a n',X' a m',-X' 

(m + n)Vri5 ntm+m , +n ,f abc f a ' b ' c 



n' .m .n,m 



16n 2 L 2 , ^ 

_L n' nm' 

9 2 L 



16n 2 L 2 , ^ 



E 



y/mnm' (m' + n') 2 
■a a V n b> a b 

a n,X a n',X' a m',-X' a m,X 

n (m! + n ) 5 m+n+n > , m i 



(17) 



(18) 



n' ,m' ,n,m 



\Jm'n'm(m + n) 
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fabc fa'b'c a' t a ] b t V /iq\ 
J J tt n',A' a n,-\ a m,\ a m'X' l iy J 



^_ g 2 L (m + n) (m! + n') 

167r2L ± n', m> n,m,k>0 Vmnm'rik 2 



5 n ,k+mfim+m> ,n+ri f f a n',\> a n,\ a m> ,\> a m,\ (20) 

32tt 2 L 2 



P* = tS^E/M<a ! <a (21) 



The self-induced inertia P 5 ~ is obtained by bringing the instantaneous inter- 
actions in Eqn. ( |i0|) into normal ordered form, Eqns. ([T7|) and (^). Before 
mass renormalization, the associated function f(m) is 

m— 1 19 Q 

/H = 8 Ep+- (MA) + 7*) + — + O (l/A) • (23) 

We choose our mass renormalization such that the logarithmic divergence is 
removed 

*t ^ o^- 1 1 3 8 4tt 2 

/(m) = 8 £p + ^ + ry (24) 

Next, we can write out terms associated with [A r , A s } 2 . Any contracted terms 
have the same operator structure as the mass term and are removed by mass 
renormalization. 

„2 T fabc fa'b'c 

p- 9 _ V J J 



167r L i k% x /klm (k + l + m) 

' a k+l+m,X'^ ' a k,X a l,X'°"m,-X (25) 



„2 T fabc fa'b'c 

p- _ 9 L, J J 

167r2L l kTL ^klm (k + l + m) 



' a m-X a 1,X' a k,X a k+l+m,X> (26) 



_2 r r fabc fa'b'c 

p- _ 9 0m+n,k+lJ J 

32tt 2 L 2 ± ^klmn 
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a fc,A'^ a LA a m,A' a n,A (27) 



„2 r r {abc fa be 

p- 9 \ - Oi +ri) k+ m J J 

10 ~~ 00^2 r 2 



-I- k.l.m.n v 



mn 



a fc,A'^ a m,-A' a LA a n,-A (28) 



„2 r r xabc fa'b'c 

p- 9 -k Ol+m,k+nJ J 

The relation between the bare mass /xo and the renormalized mass /x in 
Eqn. (|D is 

In the subsequent numerical work, we choose 



4ttL 2 



(31) 



which sets the mass scale. For the three and four particle truncations where 
we need to specify the number of colors we set N = 3. Also, we define 
M 2 as the eigenvalue of the (invariant mass) 2 operator P^P^ = 2P + P~. 
If we define the integer valued "harmonic resolution" K = LP + /(2tt), the 
continuum limit is achieved by taking lim.fr — > oo for fixed P + . States are 
classified by the conserved quantities: helicity J Z1 particle number mod 2, 
charge conjugation C, and parity V. We will often use the notation J z vc . 
Under charge conjugation, gauge fields transform as (A-)ij — * — {Ar)ji an d 
under parity, they transform as A r (x + ,x~) — > — A r (x~,x + ) (actually, P + 
eigenstates are eigenstates of parity times a longitudinal boost). Note that 
we use the definition of parity appropriate for 3+1 dimensions. Our method 
for measuring parity will be explained in more detail in Section 01. 



3 Two particles 

The low energy states in this type of model have, to a good approximation, 
definite particle number [[J. Thus, it is useful to study the spectrum by 
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Figure 1: Two particle spectrum vs fi 2 for J z — (dashed curves) and J z = 2 
(solid curves) using the polynomial wave function basis calculation discussed 
at the end of Sec. §. For J z = 0, the states alternate between even and odd 
parity. 

applying a truncation in particle number. In this section, we apply a two 
particle truncation and study the bound state equation using DLCQ and a 
polynomial wave function basis in the continuum limit. 
Two particle eigenstates have the general form 

E^SV^A^IO) (32) 

3=1 

where the total helicity is J z = A + A'. These states are even under charge 
conjugation C = 1. For J z = 2, the wavefunction must be even, <pj = (fix-j, 
since the particles are identical ^J^ c = 2 ++ ), while they can be even or 
odd for J z = [jj* = ++ or 0~ + ). The spectrum as a function of mass 
is shown in Fig. [I]. The states are roughly evenly spaced due to the linear 
potential generated by the instantaneous interactions P{ + • • • + P 4 ~. In 
addition, for n 2 = 0, it is evident that the ground state is massless in the 
continuum limit. The corresponding wavefunctions at fi 2 = are roughly 
sines and cosines as shown in Figs. |2| and |3|. 

Taking the continuum limit K — > oo, one can convert the bound state 
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Figure 2: Lowest eigenfunctions <f)j vs momentum fraction x = j/K for 
J z = 0, /i 2 = 0, and K = 100. Wavefunctions are either even or odd, 
depending on the parity of the state. 
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Figure 3: Lowest eigenfunctions <f)j vs momentum fraction x = j/K for 
J z = 2, fi 2 = 0, and K = 100. The shape of the lowest wavefunction is a 
finite K artifact; it becomes constant in the continuum limit K — > oo. 
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equation into an integral equation 

M 2 <f)(x) = fi 2 (p(x) (- + ' 
\x 



g 2 N f dy (x + y)(2-x-y) 

ny) 



4nL 2 ± 



(x - y) 2 ^jx{l - x)y{l - y) 



t£^t / = % = -M y) + 0(1/ K) . (33) 

where x,y G (0,1) and we identify <j)(i/K) with fa in the discrete case. 
For J z = 2, the coefficient of the last term of ([33|) is negative and (f>(x) = 
0(1 — x). For J z = 0, the coefficient of the last term is positive and <p(x) = 
±0(1 — x). The renormalization prescription we introduced for the discrete 
case Eqn. ( j30|) is equivalent to using a principal-value prescription for the 
pole in the first integral of Eqn. (|3^). One finds that <p(x) = 1, M 2 = 
is an exact solution for J z — 2 and fi 2 = 0, confirming our numerical 
results. Below this mass, the spectrum becomes tachyonic independent of 
the [A r ,y4 s ] 2 interaction. Also, we can repeat 't Hooft's end-point analysis 
[0, for this state. Assume that the wavefunction has power law behavior 
4>(x) oc x 13 near the endpoint x — 0. To leading order in x, the bound state 
equation becomes 

2 0_i 9 2 N /-i dy (x + y)y p - 1 / 2 



27rLj_ Jo (x — yY^Jx 

which yields the relation 

jj, 2 = i^ vr ^tan(7r ) 5) , 0</3<l/2. (35) 

7rL]_ 

Using this expression for the endpoint behavior, one can also show that M 2 
depends linearly on /i near /i 2 = 0: 



lim M = Ana, 

u2^0 



9 2N fort 
\ 47Zi • (36) 

Next, let us examine the convergence to the continuum limit for DLCQ. 
If we plot the spectrum as a function of harmonic resolution, we seem to find 
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Figure 4: Two particle spectrum for helicity 2, fi 2 = 0, as a function of 
harmonic resolution. According to this plot, convergence to the continuum 
limit seems to be quite good. 
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1/ln K 

Figure 5: Ground state eigenvalue for helicity 2, fx 2 = as a func- 
tion of harmonic resolution together with a numerical fit to 78.08/ In K — 

150.52/(lnK) 2 + 180.21/(ln_ft') 3 + In fact, convergence to the continuum 

limit is quite poor due to our having neglected the zero mode. 
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Figure 6: Lowest eigenfunction vs /i 2 for J z — 2 and K = 102. The wave- 
function is constant at jj 2 = and becomes a delta function at jj 2 = oo. 

quick convergence as shown in Fig. |4j. However, if we examine the ground 
state more closely in Fig. ||, we see that its eigenvalue in fact converges quite 
slowly. The same is true for the eigenfunctions. Let us take a variational 
ansatz for the wavefunction: <pj = jP(K — j) 13 . For J z = 2, /i 2 = 0, we 
find that this ansatz is a very good approximation of the ground state wave- 
function for all K (M 2 eigenvalues differ by only about 0.002) and that the 
wavefunction is constant (/3 — > 0) in the K — > oo limit: 

(3 w 1.71/ \nK - 2.95/(lnf0 2 + 5. 49/ (In K) 3 + . . . . (37) 

This explains the shape of the lowest eigenfunction in Figs. |3] and || 

This slow convergence is due to our having treated the zero mode sector of 
the theory poorly: we have simply removed the zero momentum mode from 
the theory. In terms of the continuum formulation, Eqn. (|33|), neglecting the 
zero mode is equivalent to introducing cutoffs in momenta. That is, 

e < x,y < 1 -e , 
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(38) 



where we identify e « 2/K. In this case the pole associated with the first 
integral is removed by adding a term of the form 

?#^> (39) 

to the right hand side of (|33|) . Using the variational ansatz <j)(x) = x^i^l — x) 13 , 
for J z = 2, fi 2 = 0, we find the same behavior that we did in the discrete 
case: 

1.838 3.399 12.611 

3 w — k o- + . . . (40 

hie (lne) 2 (hie) 3 V ; 

/iln 85.15 197.96 518.98 

me (me) (hie) 

(42) 

This slow convergence problem gradually disappears as \j? is increased since 
the amplitude of the wavefunctions becomes small near x = and 1; see 
Fig.|. 

In light of the problems associated with DLCQ, we introduce a polynomial 
wave function basis 0] to calculate the spectrum using the continuum limit 



of the bound state equation (|33"D. We assume the wavefunction is of the form 
<p{x) = x@(l — x) l3 P n (x), where P n (x) is a polynomial of order n and 3 is 
given by Eqn. (^). This method is extremely accurate. For example, with 
only five polynomials in the basis we obtain a relative error of between 10 -4 
(yU 2 w 0) and 10~ 3 (/i 2 ~ 10) for the ground state eigenvalue. 



4 Three particles 

For the three and four particle truncations, we use only the DLCQ approach 
since it would be difficult to implement a wave function basis. Unfortunately, 
the slow convergence problem that we saw in the two particle truncation will 
affect the three and four particle spectra as well. In fact, the results we 
present for finite K and \j? = better reflect the behavior of the model for 
small but nonzero /x 2 in the continuum limit. That is, the dominant finite K 
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Figure 7: Three particle spectrum for various values of J z vc where /i 2 = 
and K = 25 {J z = 1) or K = 35 {J z = 3). States whose parity could not 
be determined are placed midway between the corresponding J z +C and J z ~ c 
columns. 
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effects can be absorbed into a finite mass renormalization. In any case, the 
numbers we present should not be considered to be precise. 

The three particle spectrum is shown in Fig. [7|. Color singlet states are 
constructed from the symmetric tensor d a b c (C = — 1 ) and the antisymmetric 



tensor f a b c (C = 1) [13|. The spectrum behaves much in the same way that 
atomic spectra do: the greater the symmetry of the wavefunction, the lower 
the energy of the state. Thus, the J z = 3 sector has lower energy than the 
J z — \ sector; C = — 1 (d a bc) states have lower energy than C = 1 (f a bc) 
states; and VC = 1 states have lower energy than VC = — 1 states. In fact, 
we find the same basic behavior in the two particle sector and, for the most 
part, in the four particle sector as well. 

Next, we discuss our method for measuring parity using the three particle 
spectrum as an example. Since the parity operator exchanges x + and x~, it is 
a dynamical operator in light-front quantization and is, in general, difficult 
to construct. In addition, parity is broken by our DLCQ cutoff and can 
only be recovered in the continuum limit. (We apply periodic boundary 
conditions in the x~ direction, but not in the x + direction.) However, in a 
theory without fermions we expect that parity will be restored automatically 
in the continuum limit. To determine the parity of a state, we first determine 
its parity for finite but large /i 2 using Hornbostel's parity transform for free 



fields fill; we then decrease the mass to the value of interest. 



Three particle eigenstates have the general form: 

(fabc Or d abc ) $K,l+m+n <film Sy ^-lll) > Jz = 1 ( 43 ) 



Km n 



and 



(fab c ord abc ) ^i+m+n^^Xi^,!*! )' t7* = 3 . (44) 



l<m<n 



In Ref. [fll] Hornbostel introduces a parity transform (times a boost) that 
is valid for free fields. In general, a particle with longitudinal momentum 
fraction Xj transforms as 

Xi -> Xi = — ± . (45) 
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0.2 0.4 
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Figure 8: Hornbostel's parity transform applied to three particle wavefunc- 
tions in Eqns. (pE3| ) and (fBj). The transform maps each momentum combi- 
nation (l,m) into another momentum combination (l,m). Dots below the 
invariant curve, Eqns. (pEB]) and (p£7D, are mapped onto dots above the curve. 
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Figure 9: Contour plots of the wavef unction 4>im vs I and m for the lowest 
J^ 3 = 3 + ~ state and various values of /i 2 . Regions of negative amplitude are 
shaded gray. Comparing this plot with Fig. ||, we see that the wavefunction 
becomes an approximate eigenstate of Hornbostel's parity operator as /i 2 is 
increased. 
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In Fig. |8| we apply this transformation to Eqns. (O) and (|4*4T). In each case 
there is a curve that is invariant under the transformation: 



m = K- l -- ] jl^K-hy J z = l (46) 

and 

m = ^l ( K ~y)- l 2> J * = 3 - ( 4? ) 

For free fields, states whose wave functions <fti m are even (odd) under this 
transform have even (odd) VC. 

Since we are studying an interacting theory, we do not expect Hornbostel's 
parity operator to work very well for small fi 2 . This is, in fact, the case. On 
the other hand, as we increase fi 2 , we see that the eigenf unctions tend to 
become approximate eigenstates of Hornbostel's parity operator. Thus, we 
determine the parity of a state by applying Hornbostel's parity operator to 
the state for large fi 2 . However, in the large mass limit fi 2 — * oo, finite K ef- 
fects become important: eigenstates become states of definite momenta and 
the expectation value of Hornbostel's parity operator vanishes. Fortunately, 
there is a mass region 200 < fi 2 < 1000 where Hornbostel's parity operator 
does work fairly well. For example, if we compare Figs. |8] and |^, we see 
that, as fi 2 increases, the state becomes more approximately an eigenstate of 
Hornbostel's parity operator. 



Finally, we examine the polarized structure functions for the lowest 1 
state in Figs. |TD| and |IT]. We see that the antiparallel gluon has a somewhat 
harder momentum distribution than the two parallel gluons. This is the 
opposite of what one would normally expect on the basis of helicity retention 
|j~4|j . Physically, this result is related to the fact that the 2 ++ glueball is 
lighter than the ++ glueball. The reason is the following: If one thinks of 
the 1 glueball as a gluon-digluon bound state, then the digluon is lighter 
when it is in a 2 ++ state than when it is in a ++ state. The J z = 2 digluon 
requires an antiparallel third gluon to yield total J z = 1, while the J z — 
digluon requires a parallel gluon. Since the J z = 2 digluon is lighter, it carries 
less momentum fraction than the J z = digluon, and hence the antiparallel 
third gluon tends to be shifted towards larger momenta than the parallel third 
gluon. Of course, this argument oversimplifies the situation somewhat since 
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Figure 10: Polarized structure functions of the lowest 1 state for K = 25 
and /i 2 = 0. 




19 



there is spin dependence in the interaction between the "digluon" and the 
"third" gluon. Furthermore, since the gluons are indistinguishable one cannot 
really separate them into a "digluon" and a "third gluon" . Nevertheless, our 
main point should still be correct: that there is a connection between the 2 ++ 
glueball being lighter than the ++ glueball and having antiparallel gluons 
that are "harder" than parallel gluons. 



5 Four particles 

When one includes Fock components with a maximum of four gluons, sev- 
eral new phenomena occur: states that were present already within the two 
particle truncation mix with states that have four particles and new states 
appear that have little or no two particle content. Although the number of 
colors N must be specified for the three particle truncation, the four particle 
sector is the first place where the value of N can have an important effect on 
the spectrum. The spectrum is shown in Fig. [1^; values of M 2 quoted later 
in this section are based on this spectrum. 

The C = — 1 sector contains states that are rather high in energy. These 
states of course have no two particle content. It is rather puzzling that the 
lowest state in each C = — 1 sector has odd parity; this violates the symmetry 
principle that we discussed in Section |[ 

Let us turn our attention to the states with significant two particle con- 



tent (marked by asterisks in Fig. 12). Surprisingly, adding extra gluons has 
only very little effect on the masses and structure functions of the lightest 
glueballs (typically only a few percent change). Similar observations had 



been made in QCDi + i flTT] , [15| . The reason for this suppression is not com- 
pletely understood but has to do with the absence of UV divergences (the 
model is superrenormalizable) and IR divergences (which are screened in 
color singlet bound states) which normally drive the production of higher 
Fock components. In combination with the energy gap (due to the creation 
of an additional J z = glueball) that has to be overcome to mix in higher 
Fock components, these effects are very suppressed for low lying glueballs. 
Of course the situation changes when we examine heavier states. In that case 
there can be significant mixing of two and four particle wavefunctions. 

An example of mixing of two and four particle wavefunctions is illustrated 
in Fig. O. At fi 2 = 0, the 2 ++ sector of the theory has two nearly degenerate 
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Figure 12: Four particle spectrum for various J Z P \ K = 12 (J7^ = 0, 2) 
or = 17 (J7 Z = 4), and /i 2 = 0. States whose parity could not be deter- 
mined are placed midway between the corresponding J z +C and J z columns. 
States with more than 25% two particle content are marked with an asterisk. 
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Figure 13: Expectation value of the number operator (JV) vs M 2 and /i 2 for 
the 2 ++ sector of four particle spectrum, K = 12. We see a classic example 
of level repulsion as the two and four particle states cross. 
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Figure 14: Binding energy of the lightest 4 ++ glueball in units of the mass of 
two 2 ++ glueballs as a function of the inverse harmonic resolution for jx 2 — 1. 



states at M 2 = 91.0 and 91.3. Both of these states have significant two 
particle content. In addition, there is a four particle state at M 2 = 85.2. 
As we increase ji 2 to 1, the two particle content of the two upper states 
decreases and the lower state becomes mainly two particles. Thus, we see 
a classic example of "level repulsion" as the two particle state crosses the 
four particle state. The strength of the interaction between the two and four 
particle states is roughly (M 2 ) 24 « 2. 

Finally, we consider the C = 1 states with predominately four particle 
content. Let us start with the lightest 4 ++ glueball (M 2 = 66) as an exam- 
ple. Numerical calculations at various values of the harmonic resolution K 
(Fig. [14]) show that its mass is slightly lower than the mass of two 2 ++ glue- 
balls; that is, decay into two widely separated 2 ++ glueballs is energetically 
forbidden and the lowest 4 ++ glueball is stable. In order to minimize numer- 
ical artifacts in the analysis of this state, we compare in the following the 
results for the 4 ++ glueball taken at K = 20 with those of the 2 ++ glueball 
at half that value (K = 10). Furthermore, since no additional gluons (e.g. a 
6 particle Fock component) were allowed, we compare with the 2 ++ glueball 
restricted to the two gluon Fock sector. 

Since its mass is extremely close to the mass of two 2 ++ glueballs, one 
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Figure 15: Gluon distributions for the 2 ++ glueball (dashed curve, K = 10, 
and /i 2 = 1) and the 4 ++ glueball (solid curve, K = 20, and /i 2 = 1). We 
have rescaled the 4 ++ gluon distribution G±{x) such that it would yield the 
same distribution as the 2 ++ glueball in the zero binding limit. 
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Figure 16: Ratio of the gluon distributions shown in Fig. nj 
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would expect that this state closely resembles a loosely bound pair of 2 ++ 
glueballs. A comparison of the gluon distributions in Figs. [15] and [16] confirms 
this effect since the distribution of the 2 ++ and the (appropriately rescaled) 
4 ++ glueball are almost identical. The most pronounced differences are: 

• a slight enhancement around 0.4 < x < 0.5 which we interpret as a 
Bose-effect: when one brings two 2 ++ glueballs together, there is an 
enhancement of the gluon distribution in the region of the maximum 
due to the symmetrization of the gluon wavefunction. 

• a tail that goes beyond x > 1 in the case of the 4 ++ glueball, which 
is kinematically forbidden for the 2 ++ glueball. This behavior is also 



reflected in the divergence of the "EMC-ratio" at x = 1 in Fig. pi. 



• a depletion for x < 0.4 and for 0.5 < x < 0.95: the momentum sum 
rule requires that the EMC-ratio drops below 1 at least once, which 
would "explain" the dip near x = 0.9. However, we have not been 
able to find a similarly simple reason to explain the drop for x < 0.4. 
See, however, the discussion of the deuteron in Ref . ||15|| which showed 
a similar behavior. 

In summary, the 4 ++ glueball has much in common with a "molecule" or an 
"atomic nucleus" in the sense that many of its properties can be explained by 
a picture of weakly bound constituents (here two 2 ++ glueballs). Numerical 
limitations did not allow us to study whether similar bound states appear 
for J z = 6,8,10,.... Intuitively one would expect this to happen since (1) 
already weak attractive forces lead to binding in 1 spatial dimension and (2) 
gluons are bosons and one would not expect a shell structure such as the one 
which occurs in QCD\ + \ [fL5[]. 

We can apply a similar analysis to the J z = 2 and J z = sectors of the 
four particle spectrum. We find that the lowest 2 ++ glueball with mainly 
four particle content (M 2 = 85) is a weakly bound state of a ++ and a 2 ++ 
glueball. Likewise the lowest ++ state with mainly four particle content 
(M 2 = 70) consists of two weakly interacting 2 ++ glueballs. However, in this 
case, the state is not bound. In fact, we can also identify the lowest ++ 
state consisting of two weakly interacting ++ glueballs at M 2 = 105. This 
state is not bound either. 

The excited states in the four particle spectrum consist mainly of pairs of 
weakly interacting glueballs which are moving at different relative momenta. 
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In fact, with possibly one exception, all of the states below M 2 < 140 have 
this interpretation. For this reason, the levels shown in Fig. ^ are well 
approximated by the DLCQ formula for two free particles of invariant mass 
M x and M 2 , 

= K (M + «L) , (48) 

and therefore the level spacing of these states is a finite K artifact. The 
physically correct level spacing would be much smaller for the bound states 
and there would be a continuum spectrum in the ++ sector. 

Because of this "continuum" of weakly interacting pairs of glueballs, it is 
somewhat difficult to distinguish any "true" four particle states. The only 
candidate that we could identify is a ++ state at M 2 = 98. If we were to 
take the N — > oo limit, this would be the lowest four particle state that 
would survive. 



6 Discussion 

How well does the tube model compare with experiment? The most reliable 
results that we have to compare to are from lattice gauge theory |I6 |. The 



most striking feature of our spectrum is the fact that the 2 ++ glueball is 
lighter than the ++ glueball. This result is surely not correct. In principle, 
the coupling associated with the [A r , A,] 2 interaction in our Hamiltonian is 
independent of the coupling associated with the instantaneous interactions 
Pi+- ■ - + P± . Thus, for instance, one might adjust these couplings indepen- 
dently to obtain a spectrum that is more physically realistic @]. However, we 
note two problems with this scheme. First, we point out the special behavior 
of the two particle spectrum at ji 2 = (namely, the constant wavefunction 
and zero eigenvalue) occurs only when the two couplings are identical. More- 
over, in order to have better fit to lattice data, for instance, we would need 
to change the sign of the [y4 r ,A s ] 2 coupling. With such a sign change, the 
spectrum may be no longer be bounded from below. This is especially a 
concern if one attempts to adjust the couplings by fitting the lowest ++ and 
2 ++ glueballs to the lattice data. 

Another concern is the slow convergence of the DLCQ calculations near 
/i 2 = 0. The simplest improvement would be to choose antiperiodic instead 
of periodic boundary conditions for the gauge fields. However, if we later 
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want to include fermions in the theory, we must choose periodic boundary 
conditions for the gauge fields. Another approach would be to solve the 
constraint equation for the zero mode of A r and substitute the result back 
into the Hamiltonian. This would produce new operators in the Hamiltonian 
and would substantially improve convergence to the continuum limit. This 
technique has already been demonstrated for 4 theory in 1+1 dimensions 
[|j"7|l . Finally, one can simply abandon DLCQ as a numerical technique and 
use a polynomial wave function basis instead; this is what we did for the two 
particle truncation. However, extending such a technique to a many particle 
calculation is relatively difficult. 

In conclusion, we have made a thorough study of the 3+1 dimensional 
tube model. In particular, we have examined the behavior of the fi 2 — > 
limit and studied the convergence of the DLCQ calculations. In addition, we 
have successfully introduced an empirical technique to measure the parity of 
a state. We have found a discrete spectrum in the two and three particle 
sectors along with a continuous spectrum in the four particle sector and have 
studied the mixing of two and four particle states. 
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